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Abstract  (250  words) 


Cavitation-induced  shock  wave,  as  might  occur  in  the  head  during  exposure  to  blast  waves,  was 
investigated  as  a  possible  damage  mechanism  for  soft  brain  tissues.  A  novel  experimental 
scheme  was  developed  to  visualize  and  control  single  bubble  cavitation  and  its  collapse,  and  the 
influence  of  this  process  on  a  nearby  tissue  surrogate  was  investigated.  The  experiment  utilized  a 
Hopkinson  pressure  bar  system  which  transmits  a  simulated  blast  pressure  wave  (with  over  and 
under  pressure  components)  to  a  fluid-filled  test  chamber  implanted  with  a  seed  gas  bubble. 
Growth  and  collapse  of  this  bubble  was  visually  recorded  during  passage  of  the  blast  wave  with  a 
high  speed  camera.  To  investigate  the  potential  for  cavitation  damage  to  a  tissue  surrogate,  local 
changes  in  strain  were  measured  in  hydrogel  slices  placed  in  various  configurations  next  to  the 
bubble.  The  strain  measurements  were  made  using  digital  image  correlation  technique  by 
monitoring  the  motion  of  material  points  on  the  tissue  surrogate.  In  one  configuration,  bubble 
contact  dynamics  resulted  in  compression  contact  (>60  ps)  followed  by  inertially-driven  tension 
(>140  ps).  In  another  configuration,  the  influence  of  local  shock  waves  emanating  from 
collapsed  bubbles  was  captured.  Large  compressive  strains  (0.25  to  0.5)  that  were  highly 
localized  (0.18  mm2)  were  measured  over  a  short  time  period  (<24  ps)  after  bubble  collapse. 
High  bubble  collapse  pressures  29  to  125  times  that  of  peak  blast  overpressure  are  predicted  to 
be  the  source  of  these  localized  shock  waves.  Consistent  with  theoretical  predictions,  these 
cavitation-based  strains  are  far  larger  than  the  strains  imposed  by  passage  of  the  simulated  blast 
wave  alone.  Finally,  the  value  of  this  experimental  platform  to  investigate  the  single  bubble 
cavitation-induced  damage  in  a  biological  tissue  is  illustrated  with  an  example  test  on  rat  brain 
slices. 

Keywords:  Cavitation;  Traumatic  brain  injury;  Agarose  hydrogel;  Shock  wave,  brain 
tissue 
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Introduction 

Blast-induced  traumatic  brain  injury  (bTBI)  has  sometimes  been  referred  to  as  the  ‘signature 
injury’  of  the  Iraq  and  Afghanistan  wars,  with  estimates  of  10  to  20%  of  returning  military 
personnel  having  suffered  from  bTBI  '.  Patients  often  present  a  complex  response  to  bTBI,  and 
persistent  symptoms  can  lead  to  long-term  disabilities.  One  key  concern  is  that  pathogenesis  of 
bTBI  and  initial  mechanisms  of  injury  are  still  unclear.  Numerous  possible  mechanisms, 
including  diffuse  axonal  injury  due  to  rapid  accelerating  -  decelerating  forces,  skull  deformation 
with  elastic  flexure,  multiple  wave  reflections,  impedance  mismatch  between  heterogeneous 
tissue  types,  and  cavitation,  have  been  proposed  2.  Brain  injury  mechanisms  based  on  cavitation 
are  rooted  on  the  concept  of  shock  loading  3  where  blast  waves  consist  of  large  overpressure 
(compressive)  and  smaller  underpressure  (tensile)  components.  Propagating  blast  waves 
traveling  through  air  impinge  on  the  head,  leading  to  a  portion  of  the  wave  being  transmitted  into 
cerebrospinal  fluid  (CSF)  and  brain  tissue.  Within  the  short  time  period  of  blast  wave 
transmission,  there  is  potential  for  cavitation  bubbles  to  form,  grow  and  collapse  within  fluid 
spaces  inside  the  head.  Cavitation  occurs  spontaneously  when  liquid  of  approximately  constant 
temperature  is  subjected  to  a  sufficiently  rapid  decrease  in  pressure.  Vapor  bubbles  form  when 
pressure  falls  below  the  saturated  vapor  pressure  as  may  occur  during  the  passage  of 
underpressure  phase  of  a  blast  wave.  CSF  carries  dissolved  gasses,  and  cavitation  can  be  induced 
within  CSF  or  other  fluid  compartments  within  the  brain  when  exposed  to  relatively  low  blast 
underpressure  4.  Nucleation  of  cavitation  bubbles  may  occur  through  heterogeneous  nucleation, 
initiating  at  junctions  of  liquid-solid  boundaries,  e.g.  particulate  surfaces.  However,  pre-existing 
gas  nuclei  can  act  as  another  source  5.  Previous  hyperbaric  decompression  and  ultrasound 
studies  suggest  that  bubbles  originate  from  microscopic  gas  nuclei  present  in  biological  tissues  6" 

9 

Spatial  displacement  associated  with  formation  and  growth  of  vapor  bubbles  with  cavitation 
could  be  sufficient  to  damage  brain  tissue.  Furthermore,  considerable  energy  is  released  with 
bubble  collapse,  and  this  generates  additional  local  shock  waves  that  can  be  of  very  high 
frequency  and  magnitude  la  u.  Bubble  collapse  pressures,  which  drive  these  local  shock  waves, 
have  been  reported  over  a  wide  range  (15  MPa  to  2.5  GPa)  u"14.  These  pressures  are  potentially 
much  greater  than  blast  overpressures.  Another  characteristic  feature  of  bubble  collapse  is  the 


3 


generation  of  a  high-speed  fluid  jet  15, 16  which  may  cause  additional  damage  as  it  impinges  on 
tissue.  Resulting  collapse  events  can  mediate  highly  localized  tissue  deformation  and  damage. 
Such  cavitation  damage  of  material  surfaces  is  well  established  since  the  original  work  of 
Rayleigh  who  first  studied  pitting  of  ship  propeller  blades  or  water  turbine  blades  17' 18  which  is 
commonly  referred  to  as  ‘cavitation  erosion’.  Thus,  collapsing  cavitation  bubbles  in  the  brain 
may  cause  local  damage  to  nearby  tissues  or  blood  vessels. 

While  limitations  of  current  imaging  techniques  make  in  vivo  cavitation  measurements 
challenging,  previous  studies  suggest  that  extensive  cavitation  caused  by  explosive  blast  might 
occur  in  the  CSF  surrounding  human  brain  19-21.  For  example,  negative  pressures  have  been 
reported  in  CSF  following  exposure  to  blast  waves  generated  by  shock  tubes  in  cadaver  and 
animal  models  22.  Other  studies  using  simplified  fluid  models  of  the  head  or  head-neck  have  also 
shown  negative  pressure  to  be  induced  at  the  coutrecoup  position  when  exposed  to  a  high 
pressure  blast  wave  directed  at  coup  3'  23,  24.  In  vitro  impact  and  blast  studies  also  show  visual 
evidence  of  cavitation  3’  25 .  TBI  impact  studies  using  transparent,  fluid-filled  head  models 
showed  visual  evidence  of  bubble  formation,  growth  and  collapse  within  CSF  spaces  24.  In  shock 
tube  blast  tests,  Goeller  et  al.  3  also  showed  cavitation  at  contrecoup  regions  that  coincided 
temporally  with  periods  of  negative  pressure.  Collapse  of  cavitation  bubbles  and  elastic  rebound 
of  the  skull  resulted  in  significant  pressure  spikes  in  CSF  fluid  regions.  Other  computational 
studies  of  blast  propagation  through  the  head  and  surrounding  CSF  spaces  predict  regions  of 
negative  pressure  in  contrecoup  regions  which  also  suggest  cavitation  may  occur  during  blasts  13, 
14, 19-21  Tissue  studies  using  ultrasound  shock  waves  also  support  cavitation  induced  damage,  e.g. 
hemorrhage  and  cellular  membrane  poration  26-28 .  In  addition,  cavitation  in  blood  vessels  has 
been  shown  to  result  in  air  emboli  that  can  reach  the  brain  causing  cerebral  infarct  29, 30. 

There  are  numerous  studies  of  bubble  dynamics  during  cavitation.  The  most  common  theoretical 
model  is  the  Rayleigh-Plesset  model  1718  Plesset  and  Chapman  were  the  first  to  calculate  liquid 
jet  formation  with  collapse  of  an  empty  spherical  bubble  due  to  the  high  surrounding  fluid 
pressure  18.  Experimental  evidence  of  cavitation  is  often  given  by  acoustic  measurements  31 
and/or  optical  imaging  n.  Pressure  amplitudes  during  cavitation  bubble  collapse  have  been 
experimentally  determined  to  be  between  15  MPa  14  to  1  GPa  12  using  pressure  measurement 
systems.  By  using  Rayleigh-Plesset  theory  with  bubble  imaging  data,  maximum  collapse 
pressures  over  a  wide  range  between  8  MPa  13  to  2.5  GPa  11  have  also  been  calculated. 
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A  fundamental  problem  in  the  study  of  cavitation  is  that  a  testing  platform  does  not  readily  exist 
to  control  bubble  growth  and  collapse.  In  this  study,  a  novel  experimental  system  capable  of 
investigating  single  bubble  cavitation  was  developed.  This  experiment  utilized  a  polymer 
Hopkinson  pressure  bar  system  to  transmit  a  simulated  blast  pressure  wave  to  a  fluid-filled  test 
chamber32.  By  controlling  the  gas  content  in  the  fluid  chamber,  either  bulk  or  single  bubble 
cavitation  tests  can  be  performed.  For  single  bubble  tests,  controlled  cavitation  was  produced  by 
introducing  a  seed  gas  bubble  which  underwent  growth  and  collapse  with  blast  wave  passage. 
Cavitation  was  then  investigated  as  a  possible  damage  mechanism  for  soft  tissues  by  quantifying 
the  dynamic  behavior  of  a  nearby  tissue  surrogate  (ballistic  gel  or  0.5%  agarose  hydrogel). 
Deformation  of  speckled  hydrogel  slices  were  recorded  with  a  high  speed  camera,  and  a  digital 
image  correlation  (DIC)  technique  was  used  to  calculate  strain  fields  during  bubble  growth  and 
collapse.  The  subsequent  response  of  the  tissue  surrogate  to  the  local  blast  wave  resulting  from 
single  bubble  collapse  was  captured  and  quantified.  Implications  for  localized  high  strain  due  to 
bubble  collapse  are  discussed.  Such  a  test  system  may  be  used  to  better  isolate  damage  and 
injury  due  to  cavitation  bubbles  and,  in  future  studies,  this  system  will  be  extended  to  biological 
tissues  such  as  brain  tissue  slices. 

Methodology 

High  strain  rate  test  system 

A  split  Hopkinson  pressure  bar  (SHPB)  is  traditionally  used  to  test  metallic  and  ceramic 
materials  at  high  strain  rates  in  the  range  of  500-5000/s  33, 34.  In  this  study,  a  polymer  SHPB 
(PSHPB)  system  32,  35-37  was  used  to  generate  lower  magnitude  stress  and  pressure  waves  that 
mimicked  shock  waves  within  the  head.  A  schematic  of  the  test  system  is  given  in  Figure  1.  The 
striker  bar  (0.0254  m  diameter  and  0.38  m  in  length)  was  launched  from  a  gas  gun  towards  the 
incident  bar  (0.0254  m  diameter  and  2.44  m  in  length).  This  impact  sent  a  stress  wave  through 
the  incident  bar,  and  upon  reaching  the  piston,  the  stress  wave  accelerated  the  piston  rapidly. 
After  the  initial  compression,  a  momentum  trap  pulled  back  the  incident  bar  in  order  to  avoid 
multiple  loading  of  the  test  cell  34.  This  loaded  the  fluid  in  the  test  chamber  and  generated  a 
pressure  profile  similar  to  a  blast  wave. 

The  90  mm  x  55  mm  x  42  mm  transparent  acrylic  test  cell  consisted  of  an  interior  rectangular 
chamber  filled  with  6  ml  of  distilled  water  to  house  the  fluid  and  gel  slices.  The  interior  chamber 
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was  connected  to  a  cylindrical  chamber  in  which  the  piston  was  fit.  The  cylindrical  chamber  was 
machined  for  a  tight  fit  with  the  piston  rod,  also  O-rings  mounted  on  the  piston  provided  a  leak- 
proof  seal.  The  pressure  inside  the  chamber  was  measured  utilizing  a  high  rate  pressure  sensor 
(113B24,  PCB  Piezotronics  Inc.,  Depew,  New  York)  with  a  sample  rate  of  100  kHz.  In  addition, 
a  high  speed  digital  camera  (Vision  Research  710,  Ametek  Company,  Wayne,  New  Jersey)  was 
aligned  above  the  test  cell  to  view  cavitation  events  with  a  frame  rate  of  80,000  frames/second. 


Figure  1.  Schematic  of  the  momentum-trap  PSHPB  experimental  set-up  used  for  cavitation  testing  of 
tissue  surrogates.  The  seed  bubble  and  hydrogel  slice  specimen  were  placed  in  the  fluid  filled  test  cell  and 
loaded  by  the  incident  bar  of  PSHPB .  The  striker  bar  was  launched  by  a  gas  gun. 


Cavitation  testing 

For  bulk  cavitation  testing,  specimen  samples  were  submerged  in  the  test  cell  and  tested  with  the 
PSHPB  system.  No  bubbles  were  visible  in  the  test  cell  before  testing.  In  test  cells  and  samples 
that  were  not  degassed  or  partially  degassed,  numerous  bubbles  were  seen  to  form,  grow  and 
collapse  within  the  test  cell. 

For  single  bubble  cavitation  tests,  the  fluid-filled  test  chamber  was  degassed  to  prevent  the 
formation  of  multiple  bubbles.  To  accomplish  this,  all  parts  of  the  test  chamber  were  submerged 
in  distilled  water  in  a  large  beaker  and  degassed  within  a  vacuum  pump  (model  150,  Precision 
scientific  Inc.,  Winchester,  Virginia)  for  10  minutes.  An  ultrasonic  bath  (model  FS20,  Fisher 
Scientific  Inc.,  Hampton,  New  Hampshire)  was  then  used  to  release  and  remove  remaining  small 
bubbles  from  all  surfaces.  A  single  seed  bubble  with  a  diameter  of  less  than  50  pm  was  picked  up 
with  a  spatula  from  outside  the  chamber  and  then  placed  along  an  inside  wall  of  the  submerged 
test  cell.  To  ensure  proper  placement,  either  a  small  piece  of  tape  was  adhered  to  the  chamber 
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wall  or  a  small  region  of  the  chamber  wall  was  coated  with  petroleum  (model  p-16,  Panef  Co., 
Milwaukee,  Wisconsin)  to  act  as  a  surface  anchor  for  the  seed  bubble,  see  Figure  2. 


Pressure  sensor 


Gel  slice 


Water 


N 

Seed  bubble 


Piston 


Figure  2.  Schematic  of  the  degassed  test  chamber  with  a  single  seed  bubble  (top  view  as  seen 
by  the  high  speed  camera). 


Gel  slice  testing 

Gel  slices  were  introduced  into  the  test  cell  to  determine  the  effect  of  cavitation  on  soft  tissue 
response.  While  initial  tests  were  conducted  using  ballistic  gel  (10%,  product  1347,  VYSE, 
Addison,  IL),  the  majority  of  tests  were  conducted  with  agarose-based  hydrogel  (0.4-0. 6%, 
Trevigel  5000  powder,  Trevigen,  Inc.  Gaithersburg,  Maryland).  Hydrogels  were  prepared  in  a 
Petri  dish  and  cut  into  approximately  11  mm  x  9  mm  x  1  mm  slices.  Agarose-based  hydrogels 
have  been  previously  used  for  brain  surrogate  studies  38,  39.  Modulus  values  for  this 
concentration  (-0.27  kPa)  are  slightly  higher  than  measured  for  brain  tissue  slices  35  ’ 3S.  To  allow 
tracking  of  deformation,  hydrogel  slices  were  speckled  with  black  ink  (Accu  stamp,  Cosco 
Industries  Inc.,  Harwood  Heights,  Illinois)  with  an  air  brush  (model  200nh,  Badger  Air-Brush 
Co.,  Franklin  Park,  Illinois). 


Two  bubble  placements  were  tested:  (i)  at  the  edge  and  (ii)  underneath  slices.  With  edge 
placement  (i),  the  seed  bubble  was  placed  near  the  edge  of  a  hydrogel  slice  at  a  distance  of 
approximately  1  mm  within  the  plane  of  the  hydrogel  slice.  This  allowed  for  quantification  of 
edge  deformation  during  bubble  growth  and  collapse.  In  the  other  test  configuration  (ii),  the  seed 
bubble  was  placed  beneath  the  hydrogel  close  to  the  center  of  the  slice.  This  placement  resulted 
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in  better  visualization  of  the  local  shock  wave  response  of  the  tissue  due  to  bubble  growth  and 
collapse. 

Digital  image  correlation  (DIC) 

Digital  image  correlation  is  a  well-established  experimental  technique  to  measure  full-field 
displacements  on  any  test  sample  which  has  been  speckled  with  dye  or  ink.  These  speckles  move 
with  the  object  upon  loading  and  by  capturing  the  time  resolved  images  of  the  deformed 
specimen,  allow  measurement  of  the  displacement  on  the  entire  specimen.  This  technique  has 
been  employed  by  Subhash  et  al.  35  to  measure  tensile  deformation  of  agarose  gel  at  various 
concentrations  as  well  as  our  previous  brain  slice  testing  32.  In  the  current  study,  high  speed 
camera  images  of  the  deforming  hydrogels  were  used  to  calculate  displacement  and  strain  of  the 
hydrogel  slices  using  the  DIC  method.  2D-DIC  was  used  to  track  the  movement  of  material 
points  between  consecutive  images.  Reference  undeformed  hydrogel  slice  images  were  divided 
into  an  evenly  spaced  grid  of  dimension  (0.05  mm  x  0.05  mm)  corresponding  to  the  pixel  size  in 
each  image.  Each  pixel  had  a  different  gray  value  intensity  due  to  speckling.  A  sum  of  squared 
differences  correlation  criterion,  zero-normalized  sum  of  square  differences  (ZNSSD)  correlation 
criterion  40,  was  employed  to  track  the  displacement  of  each  pixel.  For  this  method,  a  29  x  29 
pixel  subset  centered  about  an  object  pixel  was  selected  within  the  reference  image,  and  its 
corresponding  location  in  the  deformed  image  at  the  next  time  point  (At=  12-20  s)  was 
determined  using  the  ZNSSD  correlation  function  to  evaluate  the  similarity  between  reference 
and  target  subsets.  Newton-Raphson  method  with  bicubic  spline  interpolation  was  employed  to 
solve  for  the  correlation  function.  Once  the  minimum  correlation  value  was  determined, 
differences  between  the  corresponding  positions  of  reference  and  deformed  subset  centers 
yielded  in-plane  displacements.  Strain  fields  were  calculated  from  displacement  fields  using  two- 
dimensional  Savitzky-Golay  digital  differentiation  41.  The  DIC  algorithm  was  implemented 
using  DIC  software  (MOIRE,  v.  0.955b,  Opticist.org)  based  on  Pan  et  al.  41 . 

Estimation  of  bubble  collapse  pressure 

The  Rayleigh-Plesset  equation  18  was  used  to  estimate  bubble  collapse  pressure  in  order  to 
quantify  the  local  shock  wave  effect.  The  simplified  model  is  based  on  an  isolated  spherical 
bubble  cavitation  in  an  infinite,  incompressible  liquid  and  does  not  account  for  surface  tension  or 
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viscous  effects.  The  dynamic  bubble  behavior  is  governed  by  the  following  differential  equation 

18 

-I  P  P 

rr+~r2  =  -  (l) 

2 

where  R  is  the  bubble  radius,  Pr  is  the  pressure  in  the  liquid  at  the  boundary  of  the  bubble,  Poo  is 
the  environmental  pressure,  and  p  is  the  density  of  surrounding  fluid.  For  each  test,  R(t)  was 
determined  from  high  speed  digital  images.  The  central  difference  method  was  used  to  calculate 
time  derivatives,  R(t)  and  R(t).  P,( t )  was  set  as  the  far-field  fluid  pressure  as  measured  by  the 
pressure  transducer.  The  time  point  of  bubble  collapse  was  also  determined  from  imaging  data  as 
the  point  at  which  the  bubble  shrank  to  its  smallest  measurable  radius. 

Results 

Single  bubble  cavitation 

Figure  3A  shows  temporal  evolution  of  a  typical  bubble  with  exposure  to  the  simulated  blast 
pressure  wave  shown  in  Figure  3B.  The  seed  bubble  was  compressed  during  the  compressive 
overpressure  phase  (a-c);  the  bubble  grew  during  the  tensile  underpressure  phase  (d-h);  the 
bubble  collapsed  (n);  then  there  was  secondary  growth  (o)  followed  by  secondary  collapse  of  the 
bubble  (p-r).  The  bubble  remained  approximately  spherical  throughout  the  test  until  collapse. 
The  time  from  when  the  bubble  started  to  grow  until  collapse  ranged  from  270  to  350  ps.  The 
corresponding  pressure  profile  measured  within  the  chamber  was  similar  to  a  blast  shock  wave 
with  an  overpressure  peak  phase  followed  immediately  by  a  negative  pressure  phase  and 
subsequent  smaller  oscillations.  Minimum  bubble  size  occurred  at  approximately  the  time  of 
maximum  overpressure,  P0  max .  Maximum  bubble  size  coincided  with  the  minimum 
underpressure,  Pu  min.  There  was  a  time  lag  (50  to  63  ps)  between  true  minimum  pressure  and 
maximum  bubble  size.  This  delay  was  attributed  to  (i)  the  time  required  for  the  pressure  wave  to 
travel  the  distance  between  the  bubble  and  sensor  and  (ii)  inertial  delay  of  bubble  growth.  Bubble 
collapse  was  noted  at  point  (n)  and  corresponded  to  disintegration  of  the  original  bubble  into 
multiple  smaller  bubbles. 
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Figure  3.  Dynamic  behavior  of  a  single  bubble  in  response  to  a  simulated  blast  wave.  (A) 
Sequence  of  images  showing  evolution  of  a  single  bubble  cavitation.  The  gas  seed  bubble  had  an 
initial  diameter  <  50  pm.  (B)  Corresponding  pressure  profile  showing  overpressure  and 
underpressure  phases.  t= 0  corresponds  to  the  start  of  the  PSHPB  experiment. 

Table  1  summarizes  pressure  and  bubble  behavior  during  PSHPB  testing.  For  this  test  system, 
P0  max  ranged  between  800  and  1200  kPa.  Pu  min  ranged  between  -50  and  -250  kPa.  The  average 
ratio  of  overpressure  to  underpressure  ( Pomwc  /  \Pumm\=  5.12  ±  0.45  SD)  was  relatively  constant. 
Overall,  increasing  magnitude  of  air  gun  pressure  increased  the  magnitude  of  Pu  min  and  resulted 
in  increasing  bubble  size,  see  Figure  4.  As  a  result  of  this  dependence,  bubble  behavior  was 
grouped  based  on  P„  into  low,  medium  and  high  groups,  see  Table  1.  Bubble  collapse  was 
consistently  observed  after  the  second  overpressure  peak  (point  1  in  Fig.  3,  85  ±  17  kPa)  as 
pressure  was  falling.  The  measurable  bubble  collapse  radius  varied  between  0.2  to  0.5  mm.  It 
should  be  noted  that  bubble  radius  at  collapse  was  difficult  to  determine  as  the  bubble  size 
approached  the  limit  of  image  resolution.  Therefore,  it  was  not  possible  to  identify  the  pressure 
at  the  exact  time  of  collapse  of  the  bubble  from  the  images.  To  overcome  this  limitation,  we  have 
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plotted  the  instantaneous  bubble  radius  R  during  its  growth  and  collapse  as  shown  in  Figure 
5(A).  A  fifth-order  polynomial  was  fit  to  this  data  from  which  one  can  calculate  the 
corresponding  R(t)  and  R(t)  functions  and  then  determine  the  boundary  pressure  Pc  in  the  fluid 
for  each  bubble  radius  based  on  Rayleigh-Plesset  relation  as  shown  in  Fig.  5(B).  Because  the 
radius  at  the  instance  of  bubble  collapse  (i.e.,  at  R= 0)  is  unknown,  we  fit  an  exponential  curve 
(Pc  =  keR )  for  this  data.  The  curve-fit  equations  for  the  growth  and  collapse  of  three  typical 
bubble  sizes  are  shown  in  Figure  5(B).  By  setting  R= 0,  we  determined  the  elusive  collapse 
pressure  Pc  at  the  time  of  bubble  collapse  i.e.,  Pc=k.  Typical  collapse  pressure  for  a  large  size 
bubble  (i.e.,  large  negative  or  tensile  pressures)  were  as  high  as  100  MPa.  For  comparison,  note 
from  Fig. 3(B)  that  a  typical  value  of  incident  blast  overpressure  in  our  set  up  was  only  around  1 
MPa,  and  the  blast  under  pressure  was  approximately  250  kPa.  Thus  the  local  pressure  at  the 
time  of  bubble  collapse  could  be  as  high  as  100  times  the  incident  over  pressure  or 
approximately  400  times  the  blast  under  pressure.  In  our  experiments,  Pc  was  determined  to  be 
29  to  125  times  higher  than  the  measured  peak  overpressure. 
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Figure  4.  Bubble  radius  variation  with  minimum  underpressure,  PUmm .  A  linear  function  was  fit 


(R2  =0.90). 
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Figure  5.  Bubble  time  history  and  calculated  collapse  pressure.  (A)  Typical  bubble  growth  and 
collapse  during  testing.  Radial  profiles  for  three  bubbles  corresponding  to  low,  medium  and  high 
underpressure  groups  are  shown.  5th  order  polynomials  were  fit  to  experimental  data  (R2  >0.95). 
(B)  Collapse  pressure  estimation  method.  Bubble  pressure  was  calculated  at  each  bubble  radius 
for  the  low,  medium  and  high  underpressure  tests  shown  in  (A).  Exponential  equations  were  fit 
and  pressure  at  r= 0  was  estimated  to  be  the  bubble  collapse  pressure,  Pc  (R2  >0.90). 
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Table  1.  Summary  of  test  parameters  and  bubble  behavior  during  single  bubble  cavitation  tests. 


Underpressure  group  (sample  size) 

Low  (n=10) 

Medium  (n=  8) 

High  (n=9) 

P u  min  (kPa) 

-86+11 

-166+39 

-315+21 

Air  gun  pressure  (kPa) 

69 

130 

207 

Pomax  (kPa) 

851+49 

943+71 

1178+22 

Over/underpressure  duration  (jlts) 

586+23/157+14 

643+31/168+22 

723+24/184+13 

Lag  time1"  (ps) 

53+2 

58+7 

60+3 

Maximum  radius  of  bubble  (mm) 

0.47+0.12 

1.08+0.15 

1.45+0.11 

Collapse  radius  (mm) 

0.21+0.11 

0.24+0.13 

0.23+0.09 

Poo  at  collapse*  (kPa) 

6.51+3.32 

3.38+7.44 

0.51+4.65 

Pc  (kPa) 

25.1xl03+4.2xl03 

68.6x103+5.3  xlO3 

123.9  xl03+10.2xl03 

P  c/  P  o  max 

-29.5 

-72.7 

-125.2 

*  Initial  phases 

f  Between  Pu  min  and  maximum  bubble  radius 
f  Chamber  pressure  measured  away  front  the  bubble 


Hydrogel  deformation:  High  pressure  cavitation  testing 

Initial  gel  studies  were  conducted  at  high  blast  pressures  using  fluid  chambers  that  were  not 
degassed.  PSHPB  testing  of  these  test  cells  resulted  in  multiple  bubble  formation  and  cavitation 
within  the  test  chamber.  Figure  6  shows  bubble  initiation  (b)  and  growth  (c  &  d)  near  a  gel  edge 
in  a  degassed  test  chamber.  These  experiments  showed  clear  evidence  of  shock  wave  generation 
with  bubble  collapse  (f)  that  resulted  in  hydrogel  indentation  (g-i)  well  after  the  bubble  has 
collapsed.  Thus,  this  test  showed  proof  of  bubble  collapse  and  local  shock  wave  deformation  at 
high  pressure;  other  tests  (next  section)  were  conducted  with  more  controlled  bubble  conditions 
and  revealed  the  effect  of  cavitation  at  lower  blast  pressures. 
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Figure  6.  Bulk  cavitation  and  tracking  of  a  single  bubble  near  a  ballistic-gel  edge.  Figure  (a)  shows  the 
incident  blast  wave  direction  traveling  from  right  to  left  (blue  arrows),  (b)  bubble  initiation  (circled), 
(c&d)  growth,  (f)  collapse,  and  the  effect  on  an  adjacent  gel  slice  (f-i).  A  local  shock  wave  following 
bubble  collapse  resulted  in  highly  localized  hydrogel  deformation  (red  arrows).  Pomax  =2800  kPa  and  P„ 
min=- 380  kPa. 


Hydrogel  deformation:  Single  bubble  at  an  edge 

Single  cavitation  bubble  tests  in  degassed  chambers  provided  more  consistent  bubble  placement 
~1  mm  from  the  side  edge  of  the  hydrogel  slice.  Single  bubble  dynamics  (growth  and  collapse) 
were  similar  to  that  described  in  the  previous  section  conducted  using  a  similar  blast  profile  (see 
Fig.  IB).  Hydrogel  deformation  with  direct  bubble  contact  was  captured  with  blast  wave 
propagation,  while  the  effect  of  local  shock  waves  from  bubble  collapse  was  not  detected. 
Measured  strains  in  the  vicinity  of  the  bubble  were  much  larger  than  those  which  occurred 
during  passage  of  the  compressive  ID  overpressure  wave  through  the  hydrogel  (<  0.002, 
detectable  strain  limit).  Figure  7  shows  typical  hydrogel  deformation  during  bubble  growth  and 
collapse.  As  the  bubble  grew,  the  hydrogel  edge  was  deformed  by  direct  contact  with  the 
expanding  bubble.  Radial  strain  ( err ),  circumferential  strain  (coo)  and  von  Mises  strain  (cvm)  are 
shown  in  Figure  7.  It  should  be  noted  that  deformation  and  strain  fields  were  not  calculated  in 
regions  where  the  bubble  overlapped  the  edge  of  the  hydrogel  (masked  regions)  since  there  was 
distortion  of  the  observed  speckle  field  due  to  the  overlapping  bubble.  As  these  regions 
experience  maximum  strain,  the  values  close  to  the  edge  could  not  be  calculated.  Overall,  large 
deformations  were  captured  and  magnitude  of  calculated  strains  decreased  with  distance  from  the 
bubble  center. 
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Figure  7.  Single  bubble  cavitation  at  the  edge  of  an  agarose  hydrogel  slice.  (Top  row)  High  speed  images 
showing  initial  seed  bubble  at  edge,  maximum  growth  and  collapse  at  selected  times;  (2nd  row)  radial 
strain,  (3rd  row)  circumferential  strain  and  (4th  row)  von  Mises  strain  at  the  same  time  points.  Regions  at 
the  edge  were  excluded  due  to  masking  by  the  impinging  bubble.  Test  chamber  was  degassed  before  gas 
seed  bubble  (<  50  pm  diameter)  was  placed  to  the  side  of  the  hydrogel  (~1  mm  gap  distance),  t- 0 
corresponds  to  the  start  of  the  PSHPB  experiment.  Simulated  blast  wave  traveled  from  right  to  left  with 
Pomax  =1260  kPa  and  Pu  1 89  kPa. 


To  show  dynamic  behavior,  radial  strain  of  a  material  point  along  the  edge  and  nearest  to  the 
bubble  is  shown  in  Figure  8.  err  increased  in  magnitude  and  became  more  compressive  during 
bubble  growth,  reaching  a  peak  value  (-0.15+0.03,  see  Table  2)  when  the  bubble  grew  to 
maximum  size  (1230+162  ps).  The  collapse  process  of  the  bubble  was  much  faster  compared  to 
the  retraction  time  of  the  hydrogel  due  to  the  inertia  associated  with  the  higher  density  of 
hydrogel.  As  the  bubble  shrank,  compressive  err  strains  of  the  hydrogel  returned  to  zero.  Due  to 
inertial  effects,  the  hydrogel  edge  continued  to  deform  past  its  original  position,  resulting  in 
tensile  err.  Another  peak  of  inertially-driven  tensile  strain  (0.12+0.02)  occurred  65+5  ps  after  the 
time  point  of  maximum  compression.  Maximum  calculated  strain  values  are  given  in  Table  2. 
For  all  tests,  the  maximum  compressive  err  was  slightly  greater  than  the  maximum  tensile  crr,  and 
the  ratio  of  strain  magnitudes  for  maximum  inertial  tensile  strain  to  maximum  compressive  strain 
was  0.87+0.05. 


15 


The  effect  of  local  shock  wave  generation  following  bubble  collapse  was  not  seen  by  changes  in 
err  for  this  bubble  configuration  (no  compressive  strains  after  bubble  collapse).  This  was  because 
either  the  edge  of  the  hydrogel  was  not  sufficiently  close  to  the  center  of  the  bubble  to  sense  the 
shock  wave,  or  the  expanding  bubble  pushed  the  entire  hydrogel  slice  away  from  the  bubble 
rather  than  inducing  local  deformation. 


Figure  8.  Inertially-driven  changes  in  radial  strain  in  the  hydrogel  during  bubble  growth  and  collapse. 
The  nearest  hydrogel  point  to  the  seed  bubble  that  was  not  covered  by  growing  bubble  was  selected.  The 
graph  shows  results  for  6  different  experiments.  Gas  seed  bubbles  (<  50  pm  diameter)  were  placed  to  the 
side  of  the  agarose  hydrogel  (~  1  mm  gap  distance).  t= 0  corresponds  to  the  start  of  the  PSHPB 
experiment. 


Hydrogel  deformation:  Single  bubble  underneath  a  tissue  slice 

To  better  capture  effects  of  the  local  shock  wave  generated  with  bubble  collapse,  a  single  seed 
bubble  was  placed  underneath  the  hydrogel  slice  to  provide  more  continuous  contact  with  the 
tissue  surrogate  during  the  simulated  blast.  Figure  9  shows  a  series  of  high  speed  images  for  a 
typical  experiment  of  bubble  growth  (b-c)  and  collapse  (d-e)  on  the  left  side  of  each  image,  and 
the  calculated  strain  field  in  the  hydrogel  at  that  moment  on  the  right  side  of  each  image.  Bubble 
growth  was  captured  due  to  transparency  of  hydrogel.  The  von  Mises  strain  field  was  determined 
to  be  radially  symmetric,  and  a  clear  circumferential  ring  pattern  evolved  with  time.  The  strain 
amplitude  at  the  center  of  the  bubble  reached  a  peak  value  with  maximum  bubble  expansion  (at 
c)  and  then  decreased  as  the  bubble  collapsed  (e).  After  the  bubble  collapsed  (f),  the  local  strain 
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again  increased  at  the  center  of  the  hydrogel.  These  large  compressive  strains  are  attributed  to  a 
localized  shock  wave  originating  from  the  collapsed  bubble  and  propagating  through  the 
hydrogel. 


Figure  9.  Typical  hydrogel  deformation  behavior  for  bubble  cavitation  underneath  a  speckled 
hydrogel  slice.  Left  half  of  each  of  the  high  speed  images  show  bubble  growth  (b  &  c)  and 
collapse  (d  &  e)  at  select  time  intervals  and  right  half  of  each  of  the  images  show  von  Mises 
strain  maps  calculated  using  DIC.  The  dark  region  corresponds  to  the  bubble  overlap  which 
causes  distortion  in  the  image  and  so  no  strains  were  calculated  in  this  region.  Large  compressive 
strains  due  to  localized  shock  waves  are  measured  at  the  bubble  center  after  collapse  (f).  Strains 
greater  than  0.2  were  measured  as  listed  in  Table  2;  the  selected  strain  range  highlighted  patterns 
over  the  time  range  selected.  t= 0  corresponds  to  the  start  of  the  PSHPB  experiment. 


Collapse  pressure(MPa) 

Figure  10.  Maximum  strain  due  to  collapsed  bubble  shock  wave.  Maximum  von  Mises  strain,  svm, 
following  bubble  collapse  is  plotted  against  the  calculated  collapse  pressure  which  drives  the  localized 
shock  wave.  Line  corresponds  to  a  linear  fit  (R2=  0.87). 
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Table  2.  Summary  of  cavitation  and  hydrogel  slice  test  parameters. 


Group  (size) 

Edge  placement 
(n=6) 

Underneath  placement 
(n=7) 

Air  gun  pressure  (kPa) 

207 

207 

P o  max  (kP cl) 

1100+113 

1088+86 

P u  min  (kP cl) 

-207+32 

-223+49 

Lag  time  (ps) 

61+7 

53+6 

Maximum  R  (mm) 

2.35+0.34 

2.54+0.63 

Pc  (kPa) 

1 16.3xl03+8.8  xlO3 

129.1  xl03+12.2  xlO3 

P  c/P o  max 

-105.7 

-118.7 

Maximum*  srr 
with  bubble  expansion 

-0.15+0.03 

-0.23+0.06 

Maximum  s„- 
inertially-driven  tension 

0.12+0.02 

— 

Maximum*  see 
with  bubble  expansion 

0.06+0.01 

0.17+0.03 

Maximum*  svm 
with  bubble  expansion 

0.17+0.05 

0.27+0.16 

Maximum  svm 
after  bubble  collapse 

— 

0.35+0.16 

Local  shock  area  t  (mm2) 

- 

0.18+0.03 

1  Maximum  strain  values  in  masked  regions  were  not  calculated 
t  Based  on  strain  threshold  of  svm  >  0. 15 


A  summary  of  strain  field  data  is  given  in  Table  2.  During  bubble  growth,  radial  strains  at  the 
bubble  edge  are  compressive.  Circumferential  strains  at  the  bubble  edge  are  tensile.  Following 
collapse,  highly  localized  compressive  deformation  resulted  in  maximum  strain  values  (svm 
ranged  from  0.21  to  0.51)  over  a  short  duration  (<24  ps)  at  the  location  of  the  bubble  center. 
These  strains  are  larger  than  those  determined  during  maximum  bubble  expansion.  These  strains 
are  also  much  greater  than  during  passage  of  the  compressive  ID  overpressure  wave  only  (below 
detectable  strain  limit).  Comparing  with  the  previous  test  bubble  configuration  (edge  placement, 
Figures  7  and  8),  these  localized  shock  strains  are  slightly  larger  but  over  a  shorter  time  period. 

For  each  test,  the  calculated  bubble  collapse  pressure  was  plotted  against  the  maximum  measured 
strain  after  bubble  collapse,  see  Figure  10.  The  local  shock  wave  emanating  from  the  collapsed 
bubble  is  driven  by  this  bubble  collapse  pressure,  and  the  corresponding  maximum  strain  was 
shown  to  clearly  increase  with  this  pressure  value. 
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Discussion 

In  this  study,  we  were  able  to  visually  capture  single  bubble  growth  and  collapse  as  well  as  the 
effect  of  this  single  bubble  cavitation  on  a  nearby  soft  tissue  surrogate.  Overpressure  and 
underpressure  were  controlled  by  gas  gun  pressure  of  a  PSHPB  system.  By  using  this  test 
system,  we  were  able  to  approximately  mimic  the  blast  wave  profile  that  would  be  transmitted  to 
intracranial  fluids.  The  high  overpressure  magnitudes  and  durations  were  within  relevant  ranges 
for  blast  injury  42, 43.  Underpressure  magnitudes  were  larger  than  those  reported  for  blasts  '• 44 
but  within  a  relevant  order  of  magnitude  for  injury.  We  were  able  to  capture  major  modes  of 
potential  cavitation  bubble-tissue  interaction  including  (i)  compression  contact  of  the  bubble 
growth  during  the  growth  phase,  (ii)  tensile  rebound  during  bubble  collapse,  and  (iii) 
compressive  strains  generated  by  local  shock  waves  after  bubble  collapse. 

Bubble  growth  and  collapse  were  monitored  in  degassed  test  chambers  (single  bubble  tests)  as 
well  as  chambers  that  were  not  degassed  (bulk  cavitation).  PSHPB  testing  in  non-degassed 
chambers  show  initiation  of  bubbles  and  bulk  cavitation  with  numerous  bubbles.  In  this  case, 
‘spontaneous’  bubbles  most  likely  nucleated  or  grew  from  microbubbles  that  were  not  initially 
visible  but  acted  as  natural  seeds.  Under  similar  test  conditions,  ‘spontaneous’  bubbles  showed 
similar  dynamic  behavior  to  bubbles  that  were  seeded,  indicating  that  seeded  single  bubble 
cavitation  tests  capture  relevant  cavitation  effects  needed  for  investigating  deformation  and 
injury  in  tissues.  However,  maximum  radii  were  smaller  for  the  ‘spontaneous’  bubbles.  Thus, 
lower  collapse  pressures  are  predicted  for  these  bubbles. 

In  single  bubble  testing,  maximum  bubble  size  directly  depended  on  the  minimum 
underpressure,  Pu  min,  applied  which  was  in  turn  directly  dependent  on  the  applied  PSHPB  air  gun 
pressure.  Increasing  magnitude  of  Pu  min  resulted  in  larger  bubble  size.  Secondary  oscillatory 
bubble  behavior  was  observed  with  coincident  oscillations  in  the  blast  pressure  profile.  Bubble 
duration  was  short  (~  300  ps),  and  collapse  was  visually  marked  by  fragmentation  of  the  single 
bubble  into  clusters  of  multiple  smaller  bubbles.  Calculated  values  for  collapse  pressure  at  the 
surface  of  the  bubble  also  showed  an  approximately  linear  dependence  with  Pu  min.  Thus, 
underpressure  determines  the  level  of  tissue-bubble  interaction  due  to  bubble  growth  and  the 
effect  of  localized  shock  waves  due  to  bubble  collapse. 
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To  appreciate  the  potential  for  cavitation  to  mediate  blast  injury,  it  is  important  to  emphasize  that 
even  though  the  underpressure  magnitude  is  significantly  lower  than  the  overpressure  magnitude 
in  the  incident  blast  wave,  the  bubble  collapse  pressure  is  potentially  much  higher  than  the 
applied  overpressure.  This  test  system  is  uniquely  suited  to  generate  parametric  conditions  and 
test  theoretical  predictions  about  cavitation  effects  on  tissue.  For  example,  in  Figure  5(B)  the 
predicted  collapse  pressures  were  up  to  the  order  of  100  times  that  of  overpressure  for  tested 
conditions,  and  the  test  system  was  able  to  capture  material  surface  deformation  and  strain  fields 
corresponding  to  local  shock  effects.  These  local  shock  waves  emanate  from  the  bubble  site  and 
decay  rapidly  over  a  short  distance  n.  Corresponding  tissue  surrogate  response  was  highly 
localized  (-0.18  mm2)  and  of  short  duration  (<24  ps).  Direct  evidence  of  local  shock  waves  was 
seen  only  when  the  tissue  surrogate  was  in  close  proximity  and  when  bubbles  were  placed 
underneath  slices.  Magnitude  of  these  localized  strains  was  found  to  be  dependent  on  seed 
bubble  distance  and  collapse  pressure.  Because  of  the  novel  nature  of  these  experiments, 
comparative  strain  values  for  soft  materials  is  lacking.  Strain  magnitude  is  also  dependent  on 
mechanical  properties  of  the  tissue  surrogate.  Localized  strains  were  large  (0.25  to  0.50)  and 
predicted  to  be  much  larger  than  the  compressive  strains  generated  at  earlier  time  points  with 
passage  of  the  simulated  blast  wave  through  the  hydrogel.  Large  strains  can  contribute  to 
localized  damage.  Such  direct  strain  measures  are  useful  since  strain  thresholds  are  often  what 
have  been  reported  for  brain  tissue  injury  46'48. 

The  test  system  can  also  be  adapted  to  study  tissue  injury  by  replacing  the  hydrogel  surrogate 
with  cultured  cells  or  tissue  slices.  Figure  11  shows  the  effect  of  cavitation  events  on  a  rat  brain 
tissue  slice  during  a  simulated  blast  with  bubble  initiation  and  growth  and  collapse.  Note  that  the 
maximum  bubble  size  occurs  at  t=1248  ps  where  the  tissue  tearing  starts  to  occur  well  after  the 
bubble  has  collapsed  at  time  t=1704  ps  when  localized  tissue  tearing  was  observed.  These 
observations  reveal  the  potential  for  collapse  events  to  deform  and  injure  soft  brain  tissues. 
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Figure  11.  Effect  of  cavitation  on  a  rat  brain  tissue  slice  which  was  speckled  and  suspended  in  artificial 
cerebrospinal  fluid.  A  simulated  blast  wave  traveled  from  right  to  left  with  Pomax  =1045  kPa  and  PUmm=  - 
203  kPa.  PSHPB  testing  showed  bubble  initiation  (b),  growth  (b-d)  and  collapse  (f).  Cavitation-based 
damage  with  bubble  collapse  was  observed  in  the  hypothalamus  region  (box).  Tissue  tearing  was 
observed  well  after  bubble  collapse  (f-h).  Prior  to  testing,  cell  viability  in  350  pm  thick  tissue  slices  were 
maintained  by  keeping  them  in  artificial  cerebrospinal  fluid  with  95%  O2  and  5%  CO2. 

Local  deformations  induced  by  direct  bubble  contact  were  also  captured.  Bubbles  placed 
adjacent  to  the  hydrogel  slices  rapidly  expanded  and  then  contracted,  producing  local 
deformations  in  the  tissue  surrogate.  This  resulted  in  compression  of  the  hydrogel  as  well  as 
tension  due  to  dynamic  rebounding  of  the  impinged  boundary.  Maximum  tensile  strains  were  on 
the  order  of  80%  of  the  maximum  compression  strain.  Comparing  with  the  other  test  bubble 
configuration  (underneath  placement),  these  compressive  and  tensile  strains  are  lower  but  over  a 
longer  time  period.  Such  a  dynamic  response  may  be  important  when  considering  soft  tissues 
which  are  more  sensitive  to  damage  in  tension  than  compression.  It  is  anticipated  that  in  counter 
coup  regions  where  cavitation  was  noted  in  previous  in  vivo  impact  studies  3- 49 ,  tissue  damage  is 
more  severe  and  this  severity  may  be  related  to  the  either  localized  shock  wave  compression  or 
the  tensile  component  of  the  strain  at  the  free  boundaries  of  the  brain  where  the  tissue  rebounds 
after  the  bubble  collapse.  The  current  experiments  provide  clear  evidence  predicting  such 
damage  to  occur  in  brain  tissue  during  IED  blasts. 
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Cavitation  characterization  relied  on  high  speed  imaging  which  was  optimized  to  provide  high 
resolution  images  at  frame  rates  required  to  capture  bubble  dynamics.  Pixel  size  resolution  was 
limited  when  measuring  initial  seed  bubble  size  and  size  of  the  collapsed  bubble.  Bubble  overlap 
was  also  an  issue  in  predicting  maximum  strains,  since  these  regions  were  masked.  DIC  strain 
mapping  also  assumed  2D  in-plane  deformation;  however,  some  out-of-plane  effects  were  noted 
when  bubbles  were  placed  underneath  the  slice.  As  a  result,  predicted  strains  were 
underestimated.  Measured  strains  are  also  smaller  than  would  be  expected  in  brain  tissues  since 
the  modulus  of  the  hydrogels  is  slightly  larger  than  for  brain  tissues  38 . 

The  novel  system  developed  for  cavitation  testing  generated  pressure  pulses  that  mimicked  a 
shock  wave  induced  by  an  explosive  blast,  i.e.  IED.  Single  bubble  cavitation  conditions  were 
created  and  tested  using  a  soft  tissue  surrogate.  High  speed  bubble  and  hydrogel  imaging  data 
were  used  to  elucidate  cavitation  effects  including  bubble  impingement  and  secondary  localized 
shock  waves  more  powerful  than  blast  overpressures.  Both  cavitation  effects  resulted  in  far 
greater  strains  than  seen  with  initial  passage  of  the  overpressure  component.  Thus  cavitation 
effects  should  be  further  investigated  in  tissues  as  a  localized  source  of  tissue  injury.  Future 
studies  will  further  focus  on  testing  effects  of  cavitation  on  brain  tissue  slices.  Computational 
studies  may  also  use  deformation  and  strain  maps  to  predict  stresses  within  the  tissue  substrate. 
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